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Abstract. We study the form of the luminosity distance as a function of redshift in 
the presence of large scale inhomogeneities, with sizes of order 10 Mpc or larger. We 
approximate the Universe through the Swiss-cheese model, with each spherical region 
described by the Lemaitre-Tolman-Bondi metric. We study the propagation of light 
beams in this background, assuming that the locations of the source and the observer 
are random. We derive the optical equations for the evolution of the beam area and shear. 
Through their integration we determine the configurations that can lead to an increase 
of the luminosity distance relative to the homogeneous cosmology. We find that this 
can be achieved if the Universe is composed of spherical void-like regions, with matter 
concentrated near their surface. For inhomogeneities consistent with the observed large 
scale structure, the relative increase of the luminosity distance is of the order of a few 
percent at redshifts near 1 , and falls short of explaining the substantial increase required 
by the supernova data. On the other hand, the effect we describe is important for the 
correct determination of the energy content of the Universe from observations. 
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1. Introduction 

The form of the luminosity distance as a function of redshift for distant supernovae 
supports the conclusion that the recent expansion of our Universe has been accelerating 
[1, 2]. In the context of homogeneous cosmology, a recent accelerating phase is also in 
agreement with the observed perturbations in the cosmic microwave background [3] . The 
simplest explanation for this phenomenon is that the cosmological constant is non-zero. 
However, the absence of acceleration at redshifts z ^ 1 implies that the required value 
of the cosmological constant is approximately 120 orders of magnitude smaller than its 
natural value in terms of the Planck scale. 

Various attempts at an alternative explanation try to build a link between 
inhomogeneities in the matter distribution in the Universe and the perceived cosmological 
acceleration. For example, there have been arguments, based on perturbative estimates, 
that the backreaction of super-horizon perturbations on the cosmological expansion is 
significant and could cause the acceleration [!]. However, the validity of this effect is 
questionable [5, 6]. 

We are interested in the importance for the problem of cosmological acceleration of 
inhomogeneities with sub-horizon characteristic scales today [7]. An intriguing fact is 
that the accelerating phase coincides with the period in which inhomogeneities in the 
matter distribution at scales 0(10) h' 1 Mpc become large, so that the Universe cannot 
be approximated as homogeneous any more at these scales. Such distances are still small 
compared to the Hubble distance ~ 3 x 10 3 h^ 1 Mpc, that sets the scale over which the 
Universe is assumed to be described by the homogeneous Friedmann-Robertson- Walker 
(FRW) metric. However, the presence of structures, such as walls or voids, with sizes up 
to 100 h~ x Mpc poses the question of whether their influence on the homogeneous solution 
is substantial. 

Because of the significant growth of such inhomogeneities at recent times, a 
perturbative treatment may not be sufficient. An exact solution of the Einstein equations, 
even for a simplified geometry, could be more useful in order to reveal an underlying 
mechanism. The Lemaitre-Tolman-Bondi (LTB) metric [8] has been employed often in 
this context [9]— [18]. It has been observed that any form of the luminosity distance as 
a function of redshift can be reproduced with this metric [9]. However, reproducing the 
supernova data requires a variation of the density or the expansion rate over distances 
O(10 3 ) h^ 1 Mpc [14]— [18] . Moreover, in order to avoid a conflict with the isotropy of the 
Cosmic Microwave Background (CMB), the location of the observer must be near the 
center of the spherical configuration described by the LTB metric. 

In this work we study the effect of inhomogeneities on the luminosity distance 
without assuming a preferred location of the observer. We model the inhomogeneities 
as spherical regions within which the geometry is described by the LTB metric. At 



The Effect of Large-Scale Inhomogeneities on the Luminosity Distance 



3 



the boundary of these regions, the LTB metric is matched with the FRW metric that 
describes the evolution in the region between the inhomogeneities. In this way, we 
model the Universe as consisting of collapsing or expanding inhomogeneous regions, while 
we preserve the notion of a common scale factor that describes the expansion of the 
homogeneous intermediate regions. Our model is similar to the standard Swiss-cheese 
model [19], with the replacement of the Schwarzschild metric with the LTB one. For this 
reason we refer to it as the LTB Swiss-cheese model. 

The transmission of a light beam in such a background can be studied through the 
Sachs optical equations [20] . These describe the expansion and shear of the beam along 
its null trajectory. Apart from the case of an FRW background, the optical equations have 
been derived by Kantowski for a Scharzschild background in his study of light propagation 
in the Swiss-cheese model [21]. We derive the equations for a general LTB background. 
We then use them in order to study light propagation in the LTB Swiss-cheese model. We 
focus on the modification of the luminosity distance as a function of redshift with respect 
to the homogeneous FRW case. We investigate which configurations lead to an increase 
of the luminosity distance, and estimate the magnitude of the effect. 

In the following section and in the appendices we derive the optical equations for 
LTB, FRW and Schwarzschild backgrounds. The expressions for the last two cases are 
known, but we summarize them for completeness. In section 3 we discuss the particular 
details of the LTB metric that we use for the description of the spherical regions of the 
LTB Swiss-cheese model. In section 4 we present numerical solutions for the cosmological 
evolution of the background metric. In section 5 we discuss light propagation and the 
modification of the luminosity distance in the inhomogeneous background. In section 6 
we estimate the magnitude of this modification. In section 7 we present our conclusions 
on the possible explanation of the form of the luminosity function in the context of 
inhomogeneous cosmology. 

2. Optical equations and luminosity distance 

The proper cross-section area A of a light beam obeys the equation 



where 6 is the expansion of the beam and A an affine parameter along the null beam 
trajectory. The symmetric and traceless shear tensor 



describes deformations of the beam. In this work we study light propagation in space- 
times with spherical spatial symmetry. In such cases, the off-diagonal elements o~2 of the 
shear tensor can be set consistently to zero. The reason is that the eigenvalues ±cr of 



dA 



26A. 



(2.1) 
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the shear tensor determine the deformation of a surface along two principal orthogonal 
axes perpendicular to the light direction. The rate of stretching in these two directions is 
given by 9 + o and 9 — a, respectively. In the case of spherical symmetry, one principal 
axis always lies on the plane determined by the null geodesic and the center of symmetry 
and the second perpendicularly to it. It is possible then to choose a reference basis that 
includes unit vectors along these principal directions. This permits us to set er 2 = 0, 
a — G\. 

The optical equations, derived by Sachs [20], determine the evolution of the beam 
expansion and shear as the light propagates in a background geometry. For completeness, 
we include a derivation of these equations in appendix A. In appendix B we derive their 
specific form in various backgrounds that are relevant for our study. 



2.1. Lemaitre-Tolman-Bondi (LTB) background 

Under the assumption of spherical symmetry, the most general metric for a pressureless, 
inhomogeneous fluid is the LTB metric [8]. It can be written in the form 

ds 2 = -dt 2 + b 2 (t, r)dr 2 + R 2 {t, r)dQ 2 , (2.3) 

where dfl 2 is the metric on a two-sphere. The function b(t, r) is given by 

where the prime denotes differentiation with respect to r, and f(r) is an arbitrary function. 
The bulk energy momentum tensor has the form 

T A B = dmg(-p(t,r), 0, 0, 0) . (2.5) 

The fluid consists of successive shells marked by r, whose local density p is time-dependent. 
The function R(t, r) describes the location of the shell marked by r at the time t. Through 
an appropriate rescaling it can be chosen to satisfy 

R(0,r)=r. (2.6) 

The Einstein equations reduce to 

M'(r) = 4nR 2 pR', (2.8) 

where the dot denotes differentiation with respect to t, and G = (IQirM 2 ) . The 
generalized mass function Ai(r) of the fluid can be chosen arbitrarily. It incorporates 
the contributions of all shells up to r. It determines the energy density through eq. (2.8). 
Because of energy conservation Ad(r) is independent of t, while p and R depend on both 
t and r. 
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Without loss of generality we consider geodesic null curves on the plane with 9 = tt/2. 
The geodesic equations are {k 1 = dx l /d\) 



dk c 



R'R' 



k 1 



RR^k 3 ) =0 
/' \ 



dX + 1 + / 

dk 1 -#,01.1 f R " 
— + 2-kk +^— 2(i + /) 

^ + 2 ^k°k 3 + 2^-k 1 k 3 



0. 



dX ' R" " ' R 
One of them can be replaced by the null condition 



R 



1 + / 

Eq. (2.11) can be integrated to obtain 

,3 _ c <t> 



(k'f + R 2 (k 3 ) 2 = 0. 



(2.9) 
(2.10) 
(2.11) 

(2.12) 

(2.13) 

(2.14) 
(2.15) 

JA dX^ 4M2 r V" ) " ■ (2 ' 16) 

It is clear from eq. (2.15) that the shear is generated by inhomogeneities, for which the 
local energy density is different from the average one. 



k 



R 2 ' 



The equations for the expansion 9 and the shear a of a beam take the form 
d9_ 

dX ~ AAP 



p ( k o) 2 -e 2 -a 2 



da 

Tx + 29 ° 



(k 3 ) 2 R 2 ( 3M(rY 
P 



AM 2 X AttR 3 J ' 
while the equation for the beam area can be written as 
1 d 2 VA 1 ,. n \2 



a 



2.2. Friedmann-Robertson-Walker (FRW) background 
The FRW metric is a special case of the LTB metric with 



R(t,r) = a{t)r 



P 



/(r) = cr 2 , c = 0,±l 



a 3 (t) 



M(r) 



Y c ' r 



The geodesic equations have the solution 



k 



o _ c t 



a(t) 

k ~ ± a(t) 2 
a 2 {t)r 2 



2U/2 



-1 r 2 



1/2 



(2.17) 
(2.18) 

(2.19) 

(2.20) 
(2.21) 
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The shear can be consistently set to zero as the r.h.s. of eq. (2.15) vanishes. 

The equation for the beam area (2.16) is most easily solved if the center of the 
coordinate system is taken at the location of the beam source. This implies that = in 
eq. (2.21). It can be easily checked that the solution of eq. (2.16) for an outgoing beam 
is 

A{X) = r 2 (A) a 2 (t(X)) Q s . (2.22) 

The constant fl s can be identified with the solid angle spanned by a certain beam when 
the light is emitted by a point-like isotropic source. We point out that, according to eq. 
(2.1), the beam expansion 9 diverges at the location of the point-like source (r = 0). 
We also obtain 



k l aJtt s + rak"Jtt s . (2.23) 



dX 

At the location of the source the second term vanishes, so that 



dy/A 



dX 



i — dv i — dt i — 

= fc^o) a(t.)y/n a = — (o) a(t s )sjn s = —(o) y/n a . (2.24) 

r=0 



If we normalize the scale factor so that a(t s ) = 1 at the time of the beam emission, we 
recover the standard expression A = r 2 Q s in flat space-time. In the following sections we 
study light propagation in more general backgrounds. We assume that the light emission 
near the source is not affected by the large scale geometry. By choosing an affine parameter 
that is locally A = t in the vicinity of the source, we can set 

d^A 



dX 

This expression, along with 



Q s . (2.25) 

A=0 



= 0, (2.26) 

A=0 V ' 

provide the initial conditions for the solution of eq. (8.17). For isotropic sources, we also 
expect the beam shear to vanish at the time of emission. 

2.3. Schwarzschild background 
The metric has the form 



ds 2 = - M - J dt 2 +[l~'f) 1 dr 2 + r 2 dn 2 , (2.27) 

where dfl 2 is the metric on a two-sphere and r s = Ai /(8nM 2 ) the Schwarzschild radius. 

Without loss of generality we consider geodesic null curves on the plane with 6 = it/ 2. 
The geodesic equations are [k l = dx l /dX) 

— + ^1-- k°k 1 = 2.28 
dX r 2 V r J 
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dk 1 
~dX 

dk 2 i q 
dX r 



and the null condition 



- ( 

Their solution is 

k° ■■ 



r 



r \1 
= 

k° ' 



r 



r 



e) 2 +r 2 (ey 



0. 



ct 1 



r 



k l = ± 



k z = °4. 



r_s_ 
r 



1/2 



The equations for the expansion 9 and the shear a of a beam take the form 

2 

3(A; 3 ) 2 r s 



de = _ d2 _ a2 

dX 
da 



+ 2da: 
dX 2 r 

while the equation for the beam area can be written as 
1 d 2 VA 



A dX 



-a 2 . 
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(2.29) 
(2.30) 

(2.31) 

(2.32) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 

(2.37) 



2.4- Luminosity distance and redshift 

In order to define the luminosity distance, we consider photons emitted within a solid 
angle Q s by an isotropic source with luminosity L. These photons are detected by an 
observer for whom the light beam has a cross-section A Q . The redshift factor is 



+ C0 fcg' 



'2.381 



because the frequencies measured at the source and at the observation point are 
proportional to the values of k° at these points. The energy flux f Q measured by the 
observer is 

L L fh 



fo 



(2.39) 



4ttD 2 l Att{1 + z) 2 A ' 
The above expression allows the determination of the luminosity distance D^asa function 
of the redshift z. The beam area can be calculated by solving eq. (8.17), with initial 
conditions given by eqs. (2.25), (2.26), while the redshift is given by eq. (2.38). 

In the case of a FRW background, substitution of eq. (2.22) in eq. (2.39) gives 
the standard expression D\ = (1 + z) 2 r 2 a 2 , where r Q corresponds to the comoving radial 
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coordinate of the observer (for the source r s = 0), and a Q to the value of the scale factor 
at the time of detection of the light signal. 

3. Spherical collapse in a LTB background 

The inhomogeneous cosmology we consider is a variation of the Swiss-cheese model [19]. 
Within a homogeneous background, described by a FRW metric, we consider spherical 
inhomogeneous regions. The metric within each region has a spherical symmetry around 
its center. In the traditional scenario, the total mass within every region is assumed 
to be concentrated at the center, so that the relevant metric is the Schwarzschild one. 
It is well known that the two metrics (FRW and Schwarzschild) can me matched on 
their common boundary without the appearance of a singular energy density there. The 
mass parameter of the Schwarzschild metric must be equal to the total energy within 
a homogeneous region of radius equal to that of the inhomogeneity. The study of light 
propagation in such a background was pioneered by Kantowski [21] and underlies the 
broad subject of gravitational lensing [23]. 

We are interested in light propagation in a background that accounts for the process 
of gravitational collapse. Our main aim is to examine how the luminosity-redshift 
relation is modified when the light crosses regions within which the deviation of the 
matter distribution from homogeneity is time-dependent. Within a FRW background we 
consider spherical regions where the space-time is described by the LTB metric (2.3). The 
assumption of spherical symmetry makes the problem tractable. In a realistic scenario 
we expect deviations from exact spherical symmetry, but the essense of the conclusions 
should remain unaffected. 

The choice of the two arbitrary functions Ai(r) and f(r) in eq. (2.3) can lead to 
different physical situations. The mass function A4(r) is related to the initial matter 
distribution. The function f(r) defines an effective curvature term in eq. (2.7). We can 
also interpret f(r) as part of the initial radial velocity of the fluid. We work in a gauge 
in which R(0,r) = r. We parametrize the initial energy density as Pi(r) = (1 + e(r)) /?o,ij 
with Pi(r) = p(0, r) and |e(r)| < 1. The initial energy density of the homogeneous 
background is p ,i = Po(0). If the size of the inhomogeneity is r , a consistent solution 
requires Am J r ° r 2 e(r)dr = 0, so that 



given for all r by the standard expression in homogeneous cosmology: Hf = po,i/(6M 2 ). 
Then, eq. (2.7) with R(0,r) = r implies that 




(3.1) 



We assume that at the initial time tj = the expansion rate Hi = R/R = R'/R' is 




(3.2) 
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The spatial curvature of the LTB geometry is 

(3) ^M) = -4^§- ( 3 - 3 ) 

For our choice of f(r) we find that at the initial time 

Overdense regions have positive spatial curvature, while underdense ones negative 
curvature. This is very similar to the initial condition considered in the model of spherical 
collapse [24]. 

When the inhomogeneity is denser near the center, we have f(r) < for r < r and 
f(r) = for r > tq. It is then clear from eq. (2.7) that, in an expanding Universe with 
increasing R, the central region will have R = at some point in its evolution and will 
stop expanding. Subsequently, it will reverse its motion and start collapsing. 

This point becomes clearer if we consider an example with e(r) = e\ > for r < r m , 
and e(r) = — ejrjj,/ (rjj — r^J for r m < r < tq. In the region r < r m we have r) = a(t)r, 
fir) = -(e lPo ^/6M 2 )r 2 and 

'a\ 2 1 (1 + ei)p ,i eiPo,i 1 



a/ 6M 2 a 3 6M 2 a 2 ' ^ 3 ' 5 ^ 

The evolution is typical of a closed homogeneous Universe with curvature proportional to 
€\. For r > r m , the effective curvature term 

remains negative, but goes to zero with r — > ro. The shells with r m < r < r$ stop 
expanding and eventually collapse, but at progressively later times. The shell with r = r 
expands forever. 

We can obtain an analytical expression for the growth of the central overdensity in 
this simple model. We define the quantity 

For < ei <C 1, the ratio of the energy density within the perturbation to the energy 
density far from it is given by the factor (1 + Q 3 . For £ ^ 1 and t ^ H^ 1 we find 

C-lf-i^t) . (3.8) 



5 V2 

At a time t 2 — e 1 3 ^ 2 H^ 1 we have (1 + C) 3 ~ 2, which means that the energy density within 
the perturbation is twice the value in the homogeneous region. The phase of gravitational 
collapse starts at a time t c ~ 1.5 1 2 . The growth of the perturbation ~ t 2 ^ 3 is in qualitative 
agreement with the behaviour predicted by the Jeans analysis for adiabatic subhorizon 
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perturbations in a matter dominated Universe. Our model can be viewed as an exact 
solution of the Einstein equations that is consistent with the behaviour expected from 
perturbation theory in the region of its applicability. It also has the same qualitative 
features as the model of spherical collapse [24]. 

This simple model has a physical interpretation for e± < as well. In this case 
the central region is an underdensity, surrounded by a shell with density larger than the 
average. The approximate analytical treatment remains the same as before. The effective 
curvature term is always positive and goes to zero for r —>■ r . The central region expands 
faster than the surrounding shell and its density drops faster than the average. On the 
other hand, the relative shell size shrinks and its energy density grows relative to the 
average density. We shall study a model with this typical behaviour in the following. 

We emphasize that our choice of the functions A4(r) and f(r) is not the only one 
possible. A general solution of the effective Friedmann equation (2.7) also depends on an 
arbitrary function ts(r) that determines the local Big Bang time. This results from the 
fact that eq. (2.7) can be integrated with respect to t for fixed r. The solution depends on 
the combination t — ts(R), with £e(r) arbitrary. In our model we have set ts(r) = 0. An 
example of a study of structure formation, exploring the full freedom of the LTB metric, 
is given in ref. [25]. 

One final remark concerns the form of the functions Pi{r) and e(r). If they are 
discontinuous (as in our simple example) the functions Ai'(r), f'(r) and R'(t,r) become 
discontinuous as well. This implies that the metric (2.3), (2.4) also becomes discontinuous. 
Of course the discontinuity of Pi(r) is only an approximation. The realistic physical 
situation involves a fast variation of the density. However, in order to avoid complications 
in the interpretation of our results, we always choose continuous functions Pi(r) and e(r) 
in the following. 



4. The model 



Our expressions simplify if we switch to dimensionless variables. We define t = tHi, 
f = t/tq, R = R/ro, where Hf = po,i/(6M 2 ) is the initial homogeneous expansion 
rate and r$ gives the size of the inhomogeneity in comoving coordinates. The evolution 
equation becomes 

R 3A4(f) /(f) 

R? 4vri? 3 R 2 ' 1 ' } 

with M = M/(p Qyi rl) and / = 6M 2 //(Po,i^o) = fl H h Hi = H^q. The dot now denotes 
a derivative with respect to t. 

The typical cosmo logical evolution in our model is displayed in figs. 1 and 2. The 
initial density Pi(r) = pi(r)/p = 1 + e(r) is constant pi = 1 + e\ in the region f < 0.25, 
constant pi = 1 + €2 in the region 0.5 < f < 0.75, and pi — 1 for f > 1. In the intervals 
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3.0- 




1 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

R(t,r)/R(t,1) 

Figure 1. The evolution of the density profile for a central overdensity surrounded by 
an underdensity . 

0.25 < f < 0.5 and 0.75 < f < 1 it interpolates linearly between the values at the 
boundaries. For fig. 1 we take e\ = 0.01, while for fig. 2 we use e\ = —0.01. The values 
of 62 are fixed by the requirement that Jq e(r)r 2 dr = 0. 

In fig. 1 we show the typical behaviour in the case of a central overdensity 
that is surrounded by an underdense region. We display the density profile at times 
i = 10, 100, 500, 1000. We normalize the energy density to that of a homogenous FRW 
background (given by pFRwif) — p(t, 1))- The initial time corresponds to the curve with 
the smallest deviation from 1, while the final to the curve with the largest deviation. We 
observe that the density contrast grows and eventually becomes 0(1). The central region 
becomes denser with time, and the underdense region emptier. The central region will 
stop expanding at a time i ~ 1500, reverse its motion and collapse towards the center. 
The same will happen to all the outer shells, but at progressively later times. The shell 
with f = 1 expands forever. One important feature is that the relative radius of the 
central region becomes smaller with time, as the expansion is slower than the average 
within this region. 

In fig. 2 we display the typical behaviour in the case of an underdensity surrounded 
by an overdense region. We display the density profile at times i = 10, 100, 500, 1000. 
As before, the deviation from 1 is smallest at the initial time and largest at the final. 
The central energy density drops relative to the homogeneous background, while the 
surrounding region becomes denser. The radius of the central underdensity grows relative 




Figure 2. The evolution of the density profile for a central underdensity surrounded by 
an overdensity. 



to the total size of the inhomogeneity, as this region expands faster than the average. 



5. Light propagation and luminosity distance 

We are interested in light propagation in the gravitational background of the configuration 
described in the previous section. We take the affine parameter A to have the dimension 
of time and we define the dimensionless variables A = HiX, k° = k°, k 1 = k 1 / Hi, 
k 3 = rok 3 . The geodesic equations (2.9)-(2.12) maintain their form, with the various 
quantities replaced by the barred ones, and the combination 1 + / replaced by H~ 2 + f. 
For geodesies going through subhorizon perturbations with Hi <^ 1 the effective curvature 
term / plays a minor roll. However, this term is always important for the evolution of 
the perturbations, as can be seen from eq. (4.1). 

We also define the dimensionless expansion 9 = 9 /Hi and shear a = a/Hi. The 
optical equations (2.14)-(2.16) take the form 



d\ 2 v ' V 4ttR 3 

1 d 2 VA 



A 



--p{k°y-a\ (5.4) 
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Figure 3. The luminosity distance Dl, relative to the case of homogeneous cosmology, 
for light beams crossing a configuration with a central overdensity and exiting at various 
angles <p . 



with p = p/po,i- The initial conditions (2.25), (2.26) become 
dVI 



dX 
A 



1 



n s = Jn s . (5.5) 

A=0 

0, (5.6) 



A=0 

with A = HfA and Q = HfVt. 

For our study we integrate the geodesic and optical equations numerically in the 
backgrounds described in the previous section. In order to get a feeling of the qualitative 
effect on the luminosity distance, we consider first a light beam that crosses one spherical 
inhomogeneous region. We assume that the light is emitted at a radius f = 1.5. Its 
trajectory is obtained by integrating eqs. (2.9), (2.11), (2.12) (in their rescaled form). Eq. 
2.10) is then automatically satisfied. The initial value of k$ is arbitrary, as it corresponds 
to the initial frequency of the emitted light. The integration of eq. (2.11) gives eq. (2.13). 
The angle of entry of the beam into the inhomogeneous region can be varied through the 
choice of c^, or its rescaled version = c^/r . For = the trajectory goes through 
the center of the spherical configuration. 

The form of the various light trajectories is not particularly illuminating. The 
important effect is the difference in light propagation within an inhomogeneous 
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Figure 4. The luminosity distance Dl, relative to the case of homogeneous cosmology, 
for light beams crossing a configuration with a central underdensity and exiting at various 
angles <fi . 



background and a homogeneous one. In order to make this comparison we consider two 
light trajectories (one in each background) that start at a common point with f = 1.5 
and eventually approach another common point with f = 1.5. The value of and 
the length of trajectories in terms of the affine parameter or the coordinate time are 
different. However, both light beams are emitted and received by comoving observers 
in the homogeneous region, with the same spatial coordinates. The difference in beam 
area and redshift is generated by the background that the beam crosses. Alternatively, 
one could choose to compare beams emitted at the same point within the homogeneous 
region, with the same c^. These would approach f = 1.5 at different values of the angular 
coordinate <fi. This choice gives very similar results to the ones we describe. 

In figs. 3, 4 we display the difference in luminosity distance and redshift between 
inhomogeneous and homogeneous backgrounds, as a function of the comoving angle O of 
the observer. The light is emitted from the point f s = 1.5, <f> s = 0, and observed at the 
point f a = 1.5, <f) . We have used Hi — 1 and taken the emission time of the beam to be 
t s = 550. The background within which the light propagates is not consistent with the 
theory of structure formation, as it includes a perturbation of horizon size with |ei| = 0.01, 
instead of the standard |ei| = 10 -5 . We choose such a background in order to make its 
qualitative effect on the beam properties more visible. In the following section, in which 
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we focus on a quantitative estimate of the effect, we consider realistic backgrounds. 

For the results displayed in fig. 3 the inhomogeneity consists of a central overdensity 
surrounded by an underdensity. For small values of <ft the light trajectory stays within 
the homogeneous region 1 < f < 1.5. As a result no difference is observed. For larger 
angles the trajectory crosses only the underdense part of the inhomogeneity. The resulting 
luminosity distance is larger than the one for a homogeneous background. At even larger 
angles the beam enters the central overdense region. The luminosity distance becomes 
smaller than the one in a homogeneous background. This behaviour is a consequence 
of the term ~ p(k ) 2 in the r.h.s. of eq. (2.16). If the effect is averaged over angles, a 
decrease of the luminosity distance is expected. 

The shear a has only a minor effect in the evolution of \fA and the determination of 
the luminosity distance. The reason is that the shear is generated by the difference of the 
local energy density p and an average density p av = .M/(47ri? 3 /3). These coincide within 
the central homogeneous region. They are different in the surrounding underdense region. 
However, the later gives an effect only for beams that stay within the overdensity a short 
time, without significant modification of the respective luminosity distance. A significant 
effect can be generated if the mass of the overdensity is concentrated near the center. 
In this case, eqs. (2.14)-(2.16) give an effect similar to that of eqs. (2.35)-(2.37) for a 
Schwarzschild geometry, with the additional influence of the expansion. It is known that 
for a Schwarzschild geometry the shear plays a minor role, unless the beam approaches 
the center at distances comparable to the Schwarzschild radius [21]. 

The relative difference in redshift is more than an order of magnitude smaller 
than the respective difference in luminosity distance. The reason is that the redshift, 
as opposed to the beam area, receives compensating contributions when entering and 
exiting the inhomogeneity. The net effect depends only on the evolution with time of the 
inhomogeneity (the reflection of the Rees-Sciama effect [26] in our formalism), and would 
be zero for a static one. We conclude that it is a good approximation to assume that 
the redshift remains unaffected by the presence of the inhomogeneity. This conclusion 
is consistent with the results of the study based on the standard Swiss-cheese model, in 
which the inhomogeneous region is modelled through the Schwarzschild geometry [21]. 

In fig. 4 we display the modification of the luminosity distance and redshift in the 
case of central underdensity surrounded by an overdensity. Again, we have used Hi = 1 
and taken the emission time of the beam to be i s = 550. We observe an effect on the 
luminosity distance that is opposite to that in the previous case of a central overdensity. 
Again, the modification of the redshift can be neglected in a good approximation. 

The comparison of figs. 3 and 4 demonstrates that the average modification 
of the luminosity distance, for a large statistical sample of light beams crossing the 
inhomogeneity at various angles, depends on the form of the configuration. The presence 
of underdense central regions, consistent with the appearance of large, approximately 
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z 

Figure 5. The luminosity distance Dl relative to the case of homogeneous cosmology, as 
a function of the redshift z, for light beams crossing diametrically several configurations 
with a central underdensity (upper curve) or overdensity (lower curve). 

spherical voids in the matter distribution of the Universe, is expected to lead to an increase 
of the average luminosity distance relative to the homogeneous case. The opposite is 
expected if the central regions are overdense. It must be emphasized that the average 
density is equal to that of the homogeneous case for both configurations displayed in figs. 
1 and 2. In this sense, these configurations do not have to appear in equal numbers in a 
realistic model of the Universe. It seems reasonable to suggest that the dominance of voids, 
as deduced from observations [22] , can be linked with increased luminosity distances. 

6. The magnitude of the effect 

The important question is whether the dominance of configurations with central 
underdensities can induce an effect that would explain the observed luminosity curves 
for distant supernovae. Our model is too basic to address this question in full detail. 
However, we can use it in order to estimate the magnitude of the expected effect in the 
real Universe. 

The biggest increase in the luminosity distance relative to a homogeneous cosmology 
is obtained if the Universe includes only void-like inhomogeneities with central 
underdensities. The effect is maximized if the light is assumed to cross these configurations 
passing through their center. The effect of the evolution of the inhomogeneities must be 
taken into account, with initial conditions consistent with the CMB spectrum (amplitude 



The Effect of Large-Scale Inhomogeneities on the Luminosity Distance 



17 



(9(10~ 5 ) at horizon crossing) and the observed matter distribution today. Within our 
model, this can be achieved by using as initial condition a configuration with e\ = —0.01 
and size relative to the horizon r Q Hi — Hi — 1/10 at some initial time tj = 0. This evolves 
similarly to fig. 2 into a configuration with a density contrast 0(1) at a time if ~ 1436, 
which we identify with the present. The size of the inhomogeneity relative to the horizon 
today is R(tf, r )Hf = R(tf, r ) = R(if, l)Hi. Our solution gives R(tf, 1) ~ 0.076, so that 
the size of the inhomogeneity becomes R(tf, tq) ~ 23 h~ l Mpc. This is of the order of the 
typical size of voids today. At times before the initial moment ii from which we follow 
the evolution, the inhomogeneity is only a small perturbation, with amplitude O(10~ 5 ) 
at horizon crossing. 

We consider light beams that pass repeatedly through the inhomogeneities we 
described above. The light is emitted at some time f s from a point with f = 1.5 within 
the homogeneous region. The initial conditions for the beam area are given by eqs. (2.25), 
(2.26). We assume that the light moves radially towards the center of the inhomogeneity 
(that has a fixed radius f = 1 in comoving coordinates), exits from the opposite side 
and finally arrives at a point with f = 1.5. Subsequently, the beam crosses the following 
inhomogeneity in a similar fashion. The initial conditions are now set by the values of VA 
and dVA/ dX at the end of the first crossing. Our assumption that the motion is radial 
produces the maximum effect. In general the crossing should take place at a random 
angle, with a smaller total increase of the luminosity distance. Of course, as time passes 
the profile of the inhomogeneities chang depicted in fig. 2. 

The total number of crossings determines the total redshift of the beam and the 
final beam area, related to the luminosity distance. We repeat the calculation for several 
starting times that result in a variation of the total number of crossings. In fig. 5 
(upper curve) we depict the resulting increase in the luminosity distance relative to the 
homogeneous case, as a function of the redshift. As we discussed in the previous section, 
the value of the redshift is essentially unaffected by the presence of the inhomogeneities. 
(The effect is smaller by more than an order of magnitude compared to the modification 
of the luminosity distance.) In fig. 5 we also depict the form of the luminosity curve if 
the inhomogeneity has a central overdensity (lower curve). In this case the luminosity 
distance is decreased relative to the homogeneous case. 

It is clear from fig. 5 that the influence of inhomogeneities on the luminosity distance 
is a small effect, at most of the order of a few %. This must be compared to the required 
increase in the luminosity distance in order to explain the supernova data, which is of 
the order of 30% at redshifts around 1. This conclusion remains valid for other profiles of 
the inhomogeneities, as long as the essential phenomenological requirements are satisfied 
(consistency with the observed large scale structure today, amplitude O(10~ 5 ) at horizon 
crossing) . 

It must be pointed out that the region of validity of our model does not extend to 
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very high redshifts. Despite the fact that the amplitude of the perturbation goes to zero 
at early times in our model, a significant increase of the luminosity distance is predicted 
even for large z. We do not expect such a phenomenon to prevail in a model that preserves 
the full spectrum of perturbations, instead of keeping only one with a characteristic scale, 
as we did. 

Obtaining an analytical estimate of the modification of the luminosity distance is not 
easy. However, an absolute upper bound can be derived from eq. (2.16). The "focusing" 
of a beam is minimized if the shear is negligible and the energy density of the cosmological 
medium is set to zero in this equation. In our model this idealized situation can be achieved 
if the central underdense regions of the inhomogeneities become totally empty after a long 
evolution, while overdense spherical thin shells develop around them. The effect of a shell, 
no matter how thin, on the beam area is not necessarily negligible. However, if we want 
to derive only an upper bound on the increase of the luminosity distance, we can set 
the energy density arbitrarily to zero in the optical equation only. Of course, the energy 
density still drives the cosmological expansion through eq. (2.7). The scale of the overall 
expansion is determined by the average energy density. In our model this is obvious from 
the expansion in the homogeneous regions outside the inhomogeneities. 

The equation that we derive in this way corresponds to the partially filled beam 
equation in FRW cosmology [27, 28], if all the matter is assumed to be concentrated in 
dense objects that are not crossed by the light trajectories. The same equation results 
from the optical equation (2.37) within the holes of the standard Swiss-cheese model 
if the shear is neglected. The luminosity distance as a function of the redshift can be 
derived analytically in this case, assuming that the background expansion is given by the 
standard Friedmann equation involving the average density [29]. When comparing it to 
the luminosity distance in a homogenous matter-dominated FRW cosmology, we find 



For redshifts near 1, this expression gives an increase in the luminosity distance by around 
10%. This is the maximum effect we could expect, if the expansion rate is governed by 
the average density through the standard Friedmann equation. 

7. Summary and conclusions 

The purpose of this study has been to determine the modification of the luminosity 
distance of astophysical objects, such as supernovae, which is generated by the appearance 
of large scale inhomogeneities. We constructed a picture different from that of the 
standard Swiss-cheese model. The inhomogeneities are still modelled as spherical regions 
within a homogeneous background, but the matter inside them is not concentrated in very 
dense objects at their center. Instead, it is continuously distributed in regions of density 




- 1 




- 1. 



(6.7) 
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below or above the average one. We assumed spherical symmetry for each one of the 
inhomogeneous configurations. The LTB metric gives the most general description within 
each spherical region if the cosmological fluid is assumed to be pressureless. Moreover, it 
provides an automatic matching with the intermediate homogeneous regions, as the FRW 
metric is a special case of the LTB one. As a result, our model automatically becomes an 
exact solution of the Einstein equations. 

One fundamental assumption in our study was that neither the light source nor the 
observer occupy a preferred position in the Universe. For this reason we assumed that 
the light is emitted and received within homogeneous regions, while it crosses one or 
many inhomogeneous ones along its path. We derived the necessary optical equations 
that describe the propagation of a light beam within a LTB background. 

The cosmological evolution within each spherical region has many realistic elements 
when compared with the Universe at scales around and above 10 h~ x Mpc. Central 
overdense regions become denser with time, with underdense spherical shells surrounding 
them. Central underdense regions turn into voids, surrounded by massive shells. The time 
evolution of these configurations is in qualitative agreement with perturbation theory and 
the spherical collapse model. 

The question we addressed is whether the luminosity distance of a light source 
may increase in such a background. We found that this is indeed possible. If 
the inhomogeneities involve a central underdensity and an overdense outer shell, 
the luminosity distance is larger on the average than if the matter was distributed 
homogeneously. The opposite happens if the inhomogeneity is denser in the middle, 
surrounded by an underdense shell. In this case the luminosity distance is smaller than in 
homogeneous cosmology. These results imply that a description of the Universe as being 
composed of large voids, with the matter being concentrated in the intermediate regions, 
could provide a basis for the explanation of the supernova data without dark energy. 

The problem with the picture emerging from our model is quantitative: The increase 
in the luminosity distance arising from void domination is too small to account for the 
supernova data. The relative increase at a redshift around 1 is of the order of a few %. 
Moroever, an absolute upper bound of around 10% can be derived. This is significantly 
smaller than the required increase, which is around 30%. 

We mention at this point that in our model the change in redshift because of the 
presence of inhomogeneities is smaller than the change in luminosity distance by more 
than an order of magnitude. As a result, we do not expect a large Rees-Sciama effect [26] 
that would be in conflict with the CMB data. It is also possible that the presence of voids 
could explain the observed large angle CMB anomalies [34]. 

The open question is whether an alternative model of large scale structure could result 
in a larger effect. There is a variety of metrics describing inhomogeneous cosmologies, so 
this cannot be excluded. In our modelling we have made some fundamental assumptions 
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that limit the possibilities. These are: 

a) The overall scale factor is determined by the average density through the standard 
Friedmann equation. 

b) The observer does not occupy a preferred position in the Universe, such as the center 
of a significant underdensity. In practice, we studied light signals that originate in a 
homogeneous region, cross a large number of inhomogeneities, and are detected by an 
observer within a homogeneous region. 

c) The inhomogeneities have a characteristic scale (9(10) h^ 1 Mpc today. Their evolution 
is consistent with the standard theory of structure formation. In particular, the density 
perturbations are O(10~ 5 ) at horizon crossing. 

The first two assumptions lead to the LTB Swiss-cheese model that we studied. The third 
one is consistent with the observed size of voids in the matter distribution of the Universe. 

If the first assumption is maintained we do not expect a significant increase of the 
luminosity distance, independently of the modelling of the inhomogeneous regions. The 
analytical upper bound of approximately 10% at redshifts around 1 would remain valid 
in all such cosmologies, as long as the overall scale factor is determined by the average 
density through the standard Friedmann equation. 

Of course, one could try to build a cosmology that deviates significantly from 
the FRW one, so that the effective Friedmann equation receives large corrections. An 
important limitation is that acceleration of the local volume expansion in the absence of 
vorticity requires p + 3p < 0, with p the local energy density and p the pressure [5]. 

It has been suggested that the average expansion rate in a given background may 
deviate from the local one, with the corrections accounting for the observed acceleration 
[6, 12, 30, 31]. However, how the averaging is reflected in the features of a light beam 
transmitted in the particular background is not clear yet. Obviously the light propagates 
within the exact local metric, for which the optical equations lead to an unambiguous 
determination of the luminosity distance. On the other hand, the notion of averaging 
becomes crucial for a complicated mass distribution, for which the exact local metric 
cannot be determined. In this case the averaging should be a method of determining the 
gross features of light propagation without considering all the details of the exact metric. 
More work is needed in order to have a clear understanding of the connection between 
averaging and beam features. For example, within our model the "backreaction" term 
appearing in the averaged Raychaudhuri equation [30] is zero. Despite that, we observe 
a modification of the luminosity distance. 

The basic features required for the averaging to induce a significant deviation from 
the standard Friedmann expansion were determined recently [32]. The average spatial 
curvature plays an important role in this respect, as it is coupled to large spatial variances 
of the local expansion rate and shear. Within our model the average spatial curvature 
is zero, so that this mechanism of "backreaction" through averaging is absent. This is 
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a consequence of the initial conditions we chose for the LTB metric. As we discussed 
in section 3, our initial conditions guarantee that the evolution of the inhomogeneities 
is consistent with the standard scenario of structure formation (our third assumption). 
Maintaining this consistency in backgrounds with large average spatial curvature (similar 
to those proposed in refs. [15, 16]) is a challenging problem that merits further study 
in the future. An exact background that preserves this consistency could help test the 
proposal of "backreaction" through averaging, by studying the exact light propagation in 
it. 

Another possibility in order to explain the supernova data is to violate our second 
assumption and place the observer in a special position in the Universe. It has been 
observed that any form of the luminosity distance as a function of redshift can be 
reproduced for an observer at the center of the LTB metric [9]. The reason is that, 
within an anisotropic geometry, the volume expansion rate results from the averaging of 
unequal expansion rates along various directions. In particular, for the LTB metric the 
radial expansion can be accelerating, while the tangential one decelerating [17]. Of course, 
the local volume expansion is always decelerating for p + 3p > 0. An observer located at 
the center of the LTB metric receives signals only along the radial direction, from which 
he may infer an accelerating expansion. 

Reproducing the supernova data in LTB models requires a variation of the density or 
the expansion rate over distances O(10 3 ) h" 1 Mpc [14]— [18] (also see [33]). This implies 
inhomogeneities at scales much larger than the ones allowed by our third assumption. 
Moreover, in order to avoid a conflict with the isotropy of the CMB, the location of 
the observer must be very close to the center of the spherical configuration described by 
the LTB metric. Typically the required density contrast is of order 0(1). A notable 
exception is discussed in ref. [15]. In the model presented there, the matter distribution 
is homogeneous, but the expansion rate varies by roughly 15% over distances O(10 3 ) h" 1 
Mpc. This results from a variation of the spatial curvature. Even though the conflict with 
the observed matter distribution is avoided in this way, the lack of consistency with the 
theory of structure formation and the requirement of a preferred location of the observer 
are the weak points of this approach. 

In all the models that reproduce the supernova data by placing the observer at the 
center of the LTB geometry the Universe is very inhomogeneous at earlier times. The 
problem of initial conditions remains open, while the consistency with the amplitude of 
the temperature fluctuations of the CMB is not obvious. In our approach we constrained 
our modelling of the inhomogeneities so that it leads to consistency with the theory 
of structure formation. Our findings indicate that the constraints do not permit the 
reproduction of the supernova data without dark energy. 

It is worth pointing out that the model we considered has zero vorticity. The presence 
of vorticity in the background geometry gives a positive contribution to the expansion of 
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the beam area [22]. Such models are difficult to study, but provide another unexplored 
possibility for the explanation of the supernova luminosity distance without dark energy. 

As a final remark we mention that the effect we studied in this work is important for 
the correct determination of the energy content of the Universe. Even in the presence of 
dark energy, the corrections to the luminosity distance arising from inhomogeneities must 
be taken into account, as they can modify significantly the deduced contributions from 
dark matter and dark energy to the total energy density [28]. 
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8. Appendix A 

For completeness, we derive in this appendix the optical equations in a form that is 
convenient for our study. We follow closely the presentation in ref. [22]. 

The equation for the deviation between two neighboring null geodesies x t (X) and 



D 2 C 



R l jkl Wkr- (8.1) 

The symbol D denotes a covariant derivative, while A is an affine parameter. A null 
geodesic is defined through the relations 

# = —, fc\,tf = 0, k% = 0. (8.2) 

a A 

We are interested in finding a pseudo-orthonormal basis of vectors along the path of 
light, in which we can describe the evolution of the cross section of a beam of light. As 
such we use k\ a null vector w l satisfying w % ki = —1, and L\,L\, which are space-like 
unit vectors, orthogonal to the light ray and to each other. We require that all these be 
parallely propagated along the path, so as to keep them orthogonal and normalized. In 
summary, we have along the path 

kk = ww = 0, LiLi = L 2 L 2 = 1, kw = — 1, (8.3) 

kLi = kL 2 = wLi = wL 2 = L X L 2 = 0. (8.4) 

We can choose the initial deviation such that a freely moving observer sees it 
perpendicular to the light ray: k£ l = 0, w£ l = 0. It is easy then to check that the 
deviation remains orthogonal to the rest of the path. As a result, £ l can be expressed as 
a linear combination of the vectors L l a 

f(A) = £ d a (A)^, (8.5) 

a=l,2 

where the scalars d a are the proper orthogonal components of the separation of the two 
light rays. Substituting this expansion into eq. (8.1) and using the orthogonality of L l a , 
we get 

^ = £ A ab d b , A ab = R i3k[ V a Wk k L\. (8.6) 

b 

The decomposition of the curvature tensor into the Ricci and Weyl tensors 

RijH = - (gikRji — guRjk — gjkRu + guRns) — - (gikgn — gugjk) R + CijM 

l o 

gives 

A ab = —Rijk^Sab + C ljkl D a Wk k L l b . (8.7) 
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The components d a obey the equation 

^ = E(^ + a ab )d 6 , (8.8) 
d\ b 

where 9 is the expansion of the beam and a ab the symmetric and traceless shear tensor 

a ab =( ai ^ ). (8.9) 
The later satisfies 

J2 a acO-cb = (el + of) $ab = O" 2 ^. (8.10) 

c 

The second derivative of eq. (8.8) implies that 

M^ +92+ ^ 6 + (^r + H- (8 ' n) 

Using the identity Y, a CijklL l a ki k h L l a = 0, we can see that the second term of (8.7) is 
traceless, as is the second term of (8.11). We have then 

^- + 9 2 + a 2 = -]-R ij k i k j (8.12) 
dX 2 

+ 29a ab = C i]kl U a Wk k L l b . (8.13) 

dX 

The components d a at neighboring positions along the path are 

d a = d a (A + 5X) = d a (A) + (eda + E^ d fc) SX. (8.14) 

The proper cross-section area A of the beam at neighboring positions along the path is 

A(X + 5X) = J ddndd 2 = J dd!dd 2 — = q^A(X). (8.15) 
Evaluating the Jacobian <9d/<9d results in 

^ = 26 A. (8.16) 
dX 

Eqs. (8.12) and (8.13) can now be written as 
1 d 2 y/A . 2 1 



+ o~ — ——Rijk l k 3 (8.17) 



sfA dX 2 2 

Aa Akl U n k j k K Ll 



dX 

In this work we study light propagation in space-times with spherical spatial 
symmetry. In such cases, the off-diagonal elements o~2 of the shear tensor can be set 
consistently to zero. The reason is that the eigenvalues ±<r of the shear tensor determine 
the deformation of a surface along two principal orthogonal axes perpendicular to the 
light direction. The rate of stretching in these two directions is given by 9 + o and 9 — a, 
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respectively. The off-diagonal elements o~2 are related to the orientation of the principal 
axes with respect to the space- like unit vectors L\, L|. In the case of spherical symmetry, 
one principal axis lies on the plane determined by the null geodesic and the center of 
symmetry and the second perpendicularly to it. It is possible then to take L\, L 2 along the 
principal directions. This leads to cr 2 = 0, which is equivalent to CijkiL\tik k L 2 = 0. We 
have checked explicitly that this condition is satisfied for all the space-times we consider 
in this work. As a result, eq. (8.13) becomes 

^ + 29a = C ijkl L\k j k k L[, (8.19) 

and eq. (8.18) 

= AC m L\Wk k L\. (8.20) 

The basic optical equations for our study are eqs. (8.12), (8.19), or equivalently eqs. 
(8.17), (8.20). 

9. Appendix B 

In the appendix we derive the form of the optical equations (8.16)-(8.20) in specific 
backgrounds. 

9.1. Lemaitre-Tolman-Bondi (LTB) background 

A pseudo-orthonormal basis of vectors along a null geodesic is 

k = (k°, k 1 , 0, k 3 ) (9.1) 

Li= (0,0,4,0) (9.2) 

Q,h) (9.3) 



R 

(hkPR+k l b k Q + hRk l b 

\ Wr ' ¥bR : 



k° + hR(2k l b + hk°R) k 1 b(h 2 R 2 + l) + 2hk°R h 2 R 2 -l\ . . 

W=\ n -, ,0, — — : (9.4) 

V 2{k 3 ) 2 R 2 2(k*fbR 2 2k 3 R 2 J V ; 

where the function h(r, t, u) in arbitrary. For these vectors to be parallely propagated 
along the light path, the function h must satisfy the differential equation 

$ = \- \2bhR(Rk° + R'k 1 ) + Rb 2 k l + k°R' . 

dX R 2 b l v / J 

The explicit form of the solution is not necessary for our purposes, as the quantities 

= (*°) 2 (9.5) 

CmhVk L l = i ^ - j (9.6) 

C l]kl L\Vk k L l 2 = (9.7) 
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do not depend on h. As remarked in the previous appendix, we can set a 2 = 0. The 
evolution of the expansion 9 and the shear a — o~\ of a beam is given by eqs. (2.14), 
(2.15) while the beam area evolves according to eq. (2.16). 



9.2. Schwarzschild background 

A pseudo-orthonormal basis of vectors along a null geodesic is 

k = (k°, k 1 , 0, k 3 ) (9.8) 
L 1= fo,0,-,o) (9.9) 



r 

fhk° k 1 ( r s \~ l hk l 

W 



2(k 3 )r 



k° (h 2 r 2 + l) + 2hrk x (l " ^) 
jfci (h 2 r 2 + l) + 2hrk° (l - -) , 0, k 3 (h 2 r 2 - l) V (9.11) 



r_s 
r 

where the function h(r, t, u) in arbitrary. For these vectors to be parallely propagated 
along the light path, the function h must satisfy the differential equation 

dh hk 1 k° ( r s \ 

+ + - 1-- =0. 
d\ r r z \ r J 

We have 

Rijtfk? = (9.12) 

C ijkl L{tfk k L[ = -^p^j (9.13) 
C m L\Wk*L l 2 = 0, (9.14) 

and we can set = 0. The evolution of the expansion 9 and the shear a — o~\ of a beam 
is given by eqs. (2.35), (2.36) while the beam area evolves according to eq. (2.37). 
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